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Abstract 

How molecular motors like Kinesin regulates the affinity to the rail protein in the process of ATP hydrolysis 
(ATP — > ADP Pi — ► ADP + Pi) remains to be uncovered. To understand the regulation mechanism, we 
investigate the structural fluctuation of KIF1A in different nucleotide states that are realized in the ATP 
hydrolysis process by molecular dynamics simulations of Go-like model. We found that a4 helix, which is a 
part of the microtubule (MT) binding site, changes its fluctuation systematically according to the nucleotide 
states. In particular, the frequency of large fluctuations of a4 strongly correlates with the affinity of KIF1A 
for microtubule. We also show how the strength of the thermal fluctuation and the interaction with the 
nucleotide affect the dynamics of microtubule binding site. These results suggest that KIF1A regulates the 
affinity to MT by changing the flexibility of a4 helix according to the nucleotide states. 
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INTRODUCTION 



The linear biological molecular motor is a nano-machine that converts the chemical energy 
produced by the ATP hydrolysis into the mechanical work such as cellular transport and cell 
divisions tlLl20- In contrast to macroscopic artificial machines, biomolecular motors work under 
a noisy environment such as the cell. In fact, the thermal fluctuation should be appreciable, since 
the free energy released by one ATP hydrolysis cycle is only ~ 20 k B T. Despite several decades 
of investigations, the detailed mechanism by which motor proteins use ATP to move along the rail 
protein (cy to- skeletal filament) is not understood. 

According to the nucleotide states that are realized in the ATP hydrolysis cycle, a linear motor 
generally has two binding modes to the rail protein: "strong" binding mode and "weak" binding 



mode [1, 



3, 



4J]. In the strong binding mode the motor attaches itself to the rail protein tightly, 
whereas in the weak binding mode the affinity of motor for rail protein gets lower and the motor 
readily detaches from the rail protein. The mechanism of switching between these two modes, 
however, is still unknown. To reveal this regulating mechanism, we focus on KIF1A (Kinesin-3) 
motor considering that a lot of experimental data have been accumulated so far both on structural 



and biochemical properties |3i 



KIF1A is a single-headed molecular motor and can move processively and unidirectionally 
along a microtubule (MT) by using ATP hydrolysis reaction [4]. The recent experiments which 
used several different nucleotide analogs (ATP analog, ADPPi analog, etc.) have revealed the 
structure of KIF1A and the equilibrium dissociation constant for MT (Kd) in each nucleotide 
states [3, ^ 5, 6, 7, 8]. It was found that during the ATP hydrolysis process KIF1A switches 
its binding modes for MT as follows: in the nucleotide free state and the ATP bound states, the 
"strong" binding mode is realized in which KIF1A attached tightly to the specific site of MT. In 
ADP-Pi state which is realized as a result of ATP hydrolysis reaction (ATP^ADPPi) on head, the 
binding modes of KIF1A for MT varies according to the relative configuration of ADP and 7-Pi. 
Just after the hydrolysis, in "the early ADP-phosphate state" where the distance between ADP 
and 7-Pi is small (about 3 X), the binding mode is still "strong". On the other hand, just before 
the phosphate (7-Pi) release from the KIF1A head, in "the late ADP-phosphate state" where the 
distance between ADP and 7-Pi is more than 10 X, KIF1 A realizes the "weak" binding modes in 
which the affinity of KIF1 A for MT becomes much lower than the "strong" binding modes. After 
the phosphate release, the ADP bound states realizes the "weak" binding mode. Here, we have the 
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following question: "Which part of KIF1 A regulates the switching of the affinity for MT according 
to the nucleotide states and what is its mechanism?" 

A recent experiment which includes the structural analysis suggested that mainly two loops 
(LI 1, L12) play an important role in controlling the binding modes of KIF1A for MT [3]. But it is 
also reported that other than two loops, there are a number of interaction sites with MT such as (35a, 
L8, (35b, (31a, (37b, a4, a5, which includes switch II region [8]. Then we have another question: 
"other than two loops (LI 1, LI 2), are there any part to regulate the strength of binding on MT?" 

For conventional kinesin which belongs to the same superfamily as KIF1A, it was reported 
that a large B-factor of the X-ray crystal structure does not correspond to a functionally important 
fluctuation Also, the normal mode analysis based on the elastic-network model has revealed 
that small fluctuations in the elastic regime are insufficient to capture the conformational change 
of kinesins, including KIF1A lioll . Furthermore, by a simulation study of a realistic lattice model 
it was reported that the conventional kinesin in ADP state (PDB ID code: 1BG2) exhibits partial 
folding/unfolding at the functionally important regions including the microtubule binding sites 
(switch II region) other than two loops (Lll, L12) [11]. These results indicate that larger-scale 
structural fluctuations beyond the elastic regime are essential to the function of Kinesin. 

The energy landscape theory of protein folding has been accepted widely in the last decade. The 
theory states that proteins have a funnel-like energy landscape toward the native structure Ill2l . ll3ll . 



he Go-like model is certainly the simplest class of model that realizes a funnel-like landscape 



[|14|1 and has successfully described the folding process of small proteins. Recently, some attempts 



have been made to simulate a larger structural change beyond the elastic regime by Go-like model 
jl5\\ . For conventional kinesin, in particular, Go-like model simulation succeeded in revealing 
the important function, such as the mechanism how the internal strain regulates the fluctuation of 

Ilia 1 1711 . To analyze the 



nucleotide binding site and how the directional stepping is controlled 
large structural fluctuations which is relevant to the mechanism of MT binding, we use one of the 
standard Go-models 111811 . 

In this article, we investigate the structural fluctuations of KIF1A in various nucleotide states 
at thermal equilibrium by means of the molecular dynamics simulation using a simple Go-like 
model [isl in order to reveal which part of KTF1 A plays an important role in switching the binding 
strength to MT. We also make simulations that includes nucleotide molecule explicitly to investi- 
gate the effect of presence of the nucleotide following the recent work on myosin by Takagi and 
Kikuchi [Q] 
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RESULTS 



Reference structure 

We focus on the structures of KIF1A in the following five intermediate states in the ATP 
hydrolysis: 1) the early ATP bound state in which the nucleotide binding pocket is still open. 
We call it as "ATP 1 state". 2) the late ATP bound state, or, pre-hydrolysis state, in which the 
nucleotide binding pocket is closed. We call it as "ATP 2 state". 3) the early ADP-phosphate state, 
or, post-hydrolysis state. We call it as "ADP-Pi 1 state". 4) the late ADP-phosphate state, or, pre Pi 
release state. We call it as "ADP-Pi 2 state. 5) "ADP state". The equilibrium dissociation constants 
K d from MT in l)-3) states are small (strong binding mode), whereas ones in 4)-5) states is large 
(weak binding mode) Q3L [4Q . 

We construct Go-like models for each of the five states. For the reference structures of the 
"native" states that is required to define Go-like model, we employ X-ray structures of KIF1 A with 
different nucleotide analogs. As "ATP 1 state" and "ATP 2 state", we use 1161 fe|] and 1VFV B20 
(PDB ID codes) structures that are realized in AMP-PCP and AMP-PNP (ATP analog) bound state, 
respectively. As "ADP-Pi 1 state" and "ADP-Pi 2 state", we use 1VFX Jl and 1VFZ Jl (PDB 
ID codes) structure which are realized in ADP-ALFx and ADP-Vi (ADP-Pi analog) bound state, 
respectively. Finally, as "ADP state", we use 1I5S QVD (PDB ID codes) structure which is realized 
in ADP nucleotide bound state. 

For each of the five intermediate states, we made simulations of about 5.0 x 10 7 ~ 1.0 x 10 8 
steps. The temperature was set lower than the folding temperature which was estimated by MD 
simulations (Tf ~ 1.4 — 1.45 for all the models). We mainly show the results at T = 1.3 in 
the following sections. Simulations were made by solving Langevin equation. The detail of the 
simulation is described in Model and Method section. 



Distance root-mean-square deviation for all C a pair (dRMSD) 



The Distance root-mean-square deviation (dRMSDjj) between i-j C a pair relative to the native- 
structure is defined as 
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FIG. 1: The simulated result of dRMSD contour plots for all C a pair, (a) the case in "ATP 1" (AMP-PCP 
bound) state which is a representative of strong binding states, (b) the case in "ADP-Pi 2" (ADP-Pi bound) 
state which is a representative of weak binding states. 



where N s is the step number of simulation and < > is the time average of the distance between 
i-j pair. 

In the five intermediate states, we investigated the dRMSD for all C a pair. We found a qualitative 
difference for the dRMSD between the strong binding modes ("ATP 1", "ATP 2", "ADP-Pi 1") and 
the weak binding modes ("ADP-Pi 2", "ADP"). Here we show two contour plots of dRMSD: Fig.Q] 
(a) is for "ATP 1 state", which is a representative of strong binding states, (b) is for "ADP-Pi 2 
state", which is a representative of the weak binding states. We omitted dRMSD for the gap regions 
(missing residues) which includes loop Lll (residues around 260-270) and L12 (residues around 
290-300) from Fig.Q] It can be seen in the figure that in the strong binding state ("ATP 1") a large 
structural fluctuations are localized at 270-290 residue (dRMSD ~ 10 X). On the other hand, in 
the weak binding state ("ADP-Pi 2"), no such large fluctuation is appreciable. The part (residue 
number: 270-290) that exhibits large fluctuations in the strong binding state corresponds to «4 
helix, which is a candidate of MT binding sites [80. This result suggests that the fluctuations of 
helix correlates to the strength of MT binding. To confirm it, next we investigate the dynamics of 
a 4 helix. 



Dynamics of helix 

We calculated two quantities: A6 a 4 is the angle between aA helix (residues 270-290) and a3 
helix (residues 170-190), which we choose as the representative of the head since the structure of 
ct3 helix is considerably stable and attached to the other part of the head tightly. The orientation 
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FIG. 2: A typical time courses of aA dynamics in strong binding state "ATP 1 state" (a) and weak binding 
state "ADP-Pi 2 state" (b). The upper figures show the time dependence of the value A6 a 4 (Eq.O which is 
the angle between aA helix and KIF1A catalytic core (a3 helix). The lower figures are the time dependence 
of the measure Q a ^ which is the fraction of native contacts between aA and the other parts of head. 

of two helices is presented in Fig. 0(a), which is a snapshot of the typical simulated structure of 
KIF1A in "ATP 1" state. A9„4 is the difference between the angle of the native structure and one 
of given conformation and is defined as 



180 

A0 a4 = 

7T 
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where the vectors r a3 , r a4 are given by r a3 = r 190 — r 170 , r a4 = r 2 go — r 2 7o respectively. 

The other quantity is the order parameter Q q4 , which is the native contact fraction between aA 
helix (residue 270-290) and the other parts formed in a given conformation. 

The typical time course of these two values are shown in Figf2](a) "ATP 1" (AMP-PCP bound) 
state and (b) "ADP-Pi 2" (ADP-Vi bound) state. It can be seen in the figure that aA helix in the 
strong binding state "ATP 1" exhibits considerably large fluctuations in the angle (|AO a 4| ~ 80°) 
intermittently, whereas in the weak binding state "ADP-Pi 2" the fluctuation of the angle is kept 
small (| AO a4 | < 15°). We call the intermittent large fluctuations of aA helix as "burst". The burst 
actually is a partial unfolding of the helix. We also find that the burst is accompanied by decrease 
of the order parameter Q a4: (~ 0.15). The native contact fraction of the whole head Qtotai hardly 
change (at most 10 % decrease) at the burst (time course is not shown). These facts suggest that 
breaking the contacts between aA and the other parts of head by the thermal fluctuation induces the 
burst at aA. The snapshot at the burst is shown in Fig. Ob), which corresponds to |A6 a 4| ~ 80° 
and Q a4 ~ 0.15). We find that a folded part of aA is located distant from the other parts of the 
head. 
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FIG. 3: The typical simulated structures of KIF1A in "ATP 1" (AMP-PCP bound) state. Fig.(a) is the 
snapshot when the native contacts between aA and the other part are almost formed {Q a A ~ 0.75). Fig.(b) 
is the snapshot when the burst (large fluctuation) occurs (on aA) (|AG a 4| ~ 80°). 

We also made similar simulations for other nucleotide states ("ATP 2", "ADP-Pi 1 ", and "ADP"), 
and observed the burst of aA helix. Then, we investigate the nucleotide state dependence of 
the burst frequency by making histogram of the contact fraction Q a4 from the simulated time course. 



Nucleotide state dependence of the frequency for burst fluctuation of aA helix 

The histograms of contact fraction Q a4 for different nucleotide states are shown in Fig. [4] (a). 
The distribution of Q ai varies according to the nucleotide state, while the peak position stays around 
~ 0.8 and ~ 0.2 irrespective to the nucleotide states. The ratio of the lower peak (Qai < 0.4) 
corresponds to the frequency of the burst in each state. The nucleotide state dependence of the 
frequency of the burst is shown in Fig|4](b). The horizontal axis represents the nucleotide states, 
which is arranged in the frequent order of the burst. The order 1) "ATP 1", 2) "ATP 2", 3) "ADP-Pi 
1", 4) "ADP", and 5) "ADP-Pi 2" roughly corresponds to the reaction sequence that occurs on the 
KIF1A head: "ATP 1", 2) "ATP 2", 3) "ADP-Pi 1", 4) "ADP-Pi 2", and 5) "ADP" Jj. 

From Fig.|4](b), we find that the frequency of the burst obtained by the MD correlates negatively 
with the experimental value of the equilibrium dissociation constants Kd |3j,|4J] both of CK6 (wild 
type KIF1A) and CK1 (a mutant KIF1 A that lacks the part of L12 which has positive charge and 
interacts with MT in the weak binding mode). Since the equilibrium dissociation constant is an 
inverse of the binding strength of KIF1A for the microtubule, the frequency of the burst of aA 
correlates positively to the binding strength for microtubule. These results indicate that aA helix, 
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FIG. 4: Fig. (a) is the histogram of the contact fraction Q a 4 in five intermediate states. Fig. (b) shows the 
nucleotide state dependence of the frequency of burst at temperature T = 1.3 (theory) and the equilibrium 
dissociation constants which is normalized by the one for ADP-Vi state K^(ADP — Vi). 

which is a part of MT binding sites, controls the affinity for microtubule by changing its flexibility 
according to the nucleotide state. 

The nucleotide state dependence of the frequency of the burst correlates positively with the 
binding strength regardless of temperature (T = 1.25, 1.3, 1.35 < Tf) as Fig. [51 while difference 
of the burst frequency between the weak binding states ("ADP-Pi 2" and "ADP") and the strong 
binding states ("ATP 1", "ATP 2", and "ADP-Pi 1") becomes more appreciable with temperature. 
For the temperature lower than T = 1.2, the burst is hardly observed regardless of the nucleotide 
states (data not shown). Thus, at the temperature not too lower than Tj the burst frequency is 
considerably different between the weak binding state and the strong binding state. 



The influence of the nucleotide molecule 



A protein like KIF1A work as a molecular machine only by binding the nucleotide molecule. 
To investigate the influence of the presence of the nucleotide to the structural fluctuation of the 
whole KIF1A protein, we also simulated the system which includes both KIF1A and nucleotide 
molecule explicitly. We employed the coarse-grained nucleotide model, which was introduced by 
Takagi and Kikuchi J19I1 . Here we show only a preliminary result. 

In Fig. [61 we can see the same tendency of the burst frequency as before; The burst frequency 
increases with the binding strength. The effect of the explicit nucleotide is to suppress the burst 
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FIG. 5: The nucleotide state dependence of the frequency of the burst of aA helix. This figure shows the 
effect of the temperature (T = 1.25, 1.3, 1.35) on the frequency of the burst. 




FIG. 6: This figure is the intermediate state dependence of the burst frequency {aA helix) at two condi- 
tions (filled square: simulation with explicit nucleotide molecule, open square: simulation without explicit 
nucleotide). This figure shows the effect that the presence of nucleotide gives the frequency of the burst. 

frequency, in other word, the flexibility of aA helix, to some extent. More detailed simulation with 
the explicit nucleotide model will be left for future work. 



DISCUSSION 



We have shown that according to the nucleotide states KIF1A changes the dynamics and the 
flexibility of the "«4 helix", which is a part of the MT binding site. In particular, we found that 
the frequency of the burst of aA correlated strongly with the equilibrium dissociation constants K d 
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that was obtained experimentally. This means that in the strong binding mode the MT binding site 
becomes more flexible than the weak binding mode. According to these results, we suggest the 
possibility that aA helix regulates the affinity to MT, although so far only two loops (Lll, L12) 
were considered to be important in switching the strength of binding Q3|, |4fl. The result also implies 
that the binding strength is regulated through the flexibility of the binding site. It is consistent 
with the recent simulation result for the conventional Kinesin by Kenzaki and Kikuchi 111 in . in 
which the ligand binding sites exhibit large structural fluctuations. We consider that the essence 
of above mechanism on regulating the binding strength for the rail protein also apply to other 
rail-motor systems, such as actomyosine and dynein, because the key motor elements that exhibits 
structural change in different nucleotide states include switch I and II commonly, the motifs that 



are structurally homo 
and exchange [21 , 



ogous to myosin and G protein regions that move upon nucleotide hydrolysis 

12211 . 

By recent experiments [4], the equilibrium dissociation constant K d of KIF1 A in the nucleotide 
free state is shown to be the same order as Kd in the ATP bound state, although the nucleotide 
free structure of KIF1A as well as conventional kinesin has not been solved. The present result 
suggests that judging from the value of Kd the MT binding site (piA) in the nucleotide free structure 
of KIF1 A should become considerably flexible. It seems to be consistent with a recent experiment 
for Kar3 (Kinesin- 14), which belongs to the same superfamily as KIF1A, by cryomicroscopy [24]. 
The experiment revealed that the switch II helix a4 "melts" in the nucleotide free state [12411 . 

In Fig. |4] (b), we see that the frequency of the burst correlates strongly with the equilibrium 
dissociation constant of CK1 (a mutant KIF1A) compared with that of CK6 (wild type KIF1A). 
Since CK1 lacks a part of L12 (K-loop), which has the positive charge and interact with MT 
strongly, we consider that the simulated result by our model in which the charge of the amino acid 
residues and the interaction with MT are not taken into account corresponds to a mutant CK1 rather 
than CK6. We also find that our result of burst frequency correlates strongly with the equilibrium 
constant Kd for conventional kinesin (KK1), which also lacks K-loop (see figure 2A in ref. J^]). 

We have also investigated the temperature dependence of the dynamics of a4 helix, and found 
that the burst frequency of the strong binding states and the weak binding states becomes appreciably 
different at temperature not too lower than the folding temperature. Thus, we suggest that the 
thermal fluctuation is important for switching the strength of binding on MT. 

In Fig. [7J we show the number of native contact pair between ct4 and core (other parts) in each 
intermediate state. This figure shows that the number of the native contact pair in the strong binding 
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FIG. 7: The intermediate state dependence of the number of native contact pair between aA and the other 
part (head core). 

states ("ATP 1", "ATP 2", and "ADPPi 1") is considerably smaller than that in the weak binding 
states ("ADP-Pi 2" and "ADP"). Thus we consider that the nucleotide state dependence of the burst 
frequency is a consequence of the nucleotide sate dependence of the number of the native contact 
pair between ct4 and core. 

We also found that the presence of nucleotide molecule suppresses the burst. Thus, the presence 
or absence of the nucleotide molecule affect the fluctuation of MT binding site which is located 
distant from the nucleotide binding site allosterically. 

MODEL & METHOD 
Go-like model 

To study large fluctuation of KIF1A in different nucleotide states we use "C a Go-like model" 
the version of Clementi et al. fll8Q. where a protein chain consists of spherical beads that represents 
C a atoms of amino acids residues connected by virtual bonds, and interactions are specified so that 
structures closer to the native structure are more stable. Explicitly, the effective energy V p , at a 
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protein conformation T is given as 



rW) = £ - &! 0) ) 2 + £ K t (o t - ef^f 

bonds angles 

+ Yl ^[(l-cos(0<-#)) 



dihedrals 
1 



+ i(l-cos3(0,-0S O) ))] 

native contact (0) (0) 



i<J-3 



non— native contact ^ 
i<5-3 ^ 

where r(°) signifies the native (reference) structure of protein. The vector = rj — ij is the 
distance between the ith and jth of C Q , where is the position of the ith C a . bi = |b, | = |r,- | is 
the virtual bond length between two adjacent C a . 9i is the ith angle between two adjacent virtual 
bonds, where cos 9i = (bj_i • bi)/(&i_i&»), and fa is the ith dihedral angle around bj. The first 
three terms of Eq.[3]provide local interactions, that is, bond length, bond angle, and dihedral angle 
interactions, respectively. On the other hand, the last two terms are interactions between non-local 
pairs that are distant along the chain. Native contact in the fourth term is defined as follows: if one 
of the nonhydrogen atoms in the zth amino acid is within a distance of 6.5 °A from any nonhydrogen 
atom in the jth amino acid, we define the pair of the ith and jth amino acids as being native 
contact. Parameters with the subscript (0) are the constants, of which values are taken from the 
corresponding variables in the native structure. Therefore, all of the terms except the last one are 
set up so that each term has the lowest energy when the conformation T coincides with the native 
structure r( ): this effect realizes the funnel-like energy landscape. 

According to nucleotide states, KIF1 A has different gap regions (missing residues) which exists 
mainly around two loop (Lll, L12) region (Lll around 260, L12 around 290) and C terminal 
(around 360) J3, 7]. The length of gap regions changes according to nucleotide state. Since in 
the gap regions there are not any information about native structure, we set the bond length bf^ 1 in 
gap region 3.8 A which is the average value of protein, and set the interaction parameters in gap 
region (K e , K^, k nc ) zero, whereas the values of the others (parameters K r , k nnc ) are hold. Because 
the residues in gap regions are assumed to fluctuate freely except restraint on the bond length, the 
above parameter sets for gap regions may be appropriate. 

If the sets for gap regions is excluded, the interactions parameter we use throughout the present 
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work K r = 100.0, K e = 20.0, Ka 



Hi 



.0, k QC 



0.25 and C = 4.0 X are the same values as 



those used in Takagi and Kikuchi [L9J]. (The cutoff length for calculating the forth term in Eq. |3] 
is also taken to be 2ry like a past work [19].) 



Dynamics 

The dynamics of protein are simulated by the underdamped Langevin equation at a constant 
temperature T (in the thermal equilibrium). 

miWi = Fi - 7iVi + (4) 

where v, is the velocity of the ith bead and a dot represents the derivative with respect to time t (thus, 
Vj = Vi). Fi and & are systematic and random forces on ith bead, respectively. The systematic force 
Fj is derived from the effective energy V p and can be defined as F* = — dV p /dr i . £j is a Gaussian 
white random forces, which satisfies (&) = and (£i{t)£j(t)} = 2jT8ij5(t — t')l, where the 
bracket denotes the ensemble average and 1 is a 3 x 3 unit matrix. Here, we note that the same unit 
is used both for energy and temperature and thus the Boltzmann constant fee = 1. For a numerical 



integration of the Langevin equation, we use an algorithm by Honeycutt and Thirumalai [20] . We 
use 7 = 0.25, mj=1.0, and the finite time step At = 0.02. 

For a given protein conformation T, we define that the native contact between % and j is formed 
if the distance is < 1.2r^\ where rf^ 1 is the distance between the ith and jth amino acids at the 
native structure I^ '. We then use a standard measure of the nativeness, Q(T), for a given protein 
conformation T, defined as the ratio of numbers of formed native contacts at T to those at the native 
structure r^ ^. 

We thank Shoji Takada and Kazuo Sasaki for valuable discussions. 
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